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Abstract 

The optical conductivity measurements give a powerful tool to investigate the nature of the 
superconducting gap for conventional and unconventional superconductors. In this article, first, 
general analyses of the optical conductivity are given stemmed from the Mattis-Bardeen formula for 
conventional BCS superconductors to unconventional anisotropic superconductors. Second, we dis- 
cuss the reflectance-transmittance (R-T) method which has been proposed to measure far-infrared 
spectroscopy. The R-T method provides us precise measurements of the frequency-dependent con- 
ductivity. Third, the optical conductivity spectra of the electron-doped cuprate superconductor 
Nd2-a;Ce2;Cu04 are investigated based on the anisotropic pairing model. It is shown that the 
behavior of optical conductivity is consistent with an anisotropic gap and is well explained by 
the formula for d-wave pairing in the far-infrared region. The optical properties of the multiband 
superconductor MgB2, in which the existence of superconductivity with relatively high- Tc (39K) 
was recently announced, is also examined to determine the symmetry of superconducting gaps. 
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I. INTRODUCTION 



The measurements of optical properties provide us important insights concerning the 
nature of charge carries, pseudogaps and superconducting g cipS, clS well as the electronic 
..ucujot a. .«a,.0THe optica, .pec..o.op; gives a view into the electronic 
structure, low-lying excitations, phonon structure, etc. The optical conductivity or the di- 
electric function indicates a response of a system of electrons to an applied field. For the 
ordinary superconductors the evidence for an energy gap has been obtained by infrared 
spectroscopy. Far above the superconducting energy gap, a bulk superconductor behaves 
like a normal metal in the optical response. The Mattis-Bardeen formula derived in the 
BCS theory consistently describes the infrared behaviors in the classical conventional super- 
conductors. After the discovery of high-temperature superconductivity, a large amount of 
works has been made to find the superconducting gap and any spectral features responsible 
to the superconducting pairing, using an infrared spectroscopy technique. 

Optical properties are discussed in the linear response theory where the induced currents 
are proportional to the external applied field. General formulas have been derived for the 
optical response. In this paper in Section II we discuss the linear response theory for the 
conductivity; we derive the Mattis-Bardeen formula for conventional superconductors and 
the formula for London superconductors. The conductivity sum rule is briefiy discussed 
here. 

In Section III we briefiy present a new method to characterize far-infrared optical prop- 
erties which we call the refiectance-transmittance method (R-T method). In this method, 
both the refiectance spectra R{uj) and the transmittance spectra T{uj) are measured, and 
then they are substituted into a set of coupled equations which describe exactly the trans- 
mittance and refiectance of thin films. The coupled equations are solved numerically by the 
Newton method to obtain the complex refractive indices n and k of thin films as functions 
of the frequency uj, which determined the optical conductivity cr(u;). Since this method 
does not need a Kramers-Kronig transformation, we are free from difficulties stemmed from 
uncertainties in the small u region in the conventional method. 

In the subsequent Sections we discuss two materials: the electron-doped oxide super- 
conductor Nd2_a;Ce2^Cu04 and the magnesium diboride MgB2 exhibiting = 39K. The 
cuprate high-Tc superconductors are regarded as a typical London superconductor satisfy- 
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ing A » ^ for the penetration depth and the coherence length. Our date obtained from the 
R-T method for Nd2_a;Cej;Cu04 clearly indicates the d-wave symmetry with nodes for the 
superconducting gap. 

MgB2 is recently discovered superconductor with a relatively high in spite of its simple 
crystal structure. The symmetry of Cooper pairs is an issue which should be clarified to 
investigate the mechanism of high Tc. The optical properties provide us information on 
superconducting gaps from which we conclude that this material is described by two order 
parameters attached to a- and 7r-bands. Besides, the two order parameters have different 
anisotropy to explain the experimental results consistently. 

II. THEORY OF OPTICAL CONDUCTIVITY 
A. Linear Response Theory 

In this section we discuss the optical properties in the linear response theory in the 
normal metal and superconductors. The famous Mattis-Bardeen formula is derived and its 
modifications to unconventional superconductors are discussed. In the Kubo theory the 
external applied field 

E\t) = -Y. ^Ait) = - E ^A^"^'^ (2-1) 

is considered as a perturbation to the non- interact ing system described by the Hamiltonian 
Hq. The total Hamiltonian is given hy H = Hq + H'. From the equation ihdp/dt = [H,p], 
a linear variation p'{t) to the density operator po = e~^^° /Zq is written as 

ih^ = [Ho,p'] + [H',Po]. (2.2) 

Then p' is given as 

p'{t) = ~ J' dt'e-'"°^'-''y^[H'{t'),po]e'''^'-''^^''. (2.3) 
For the electrical conductivity, the external fields are given by 

a^^e^Xiu,, X^ = E^ {p^x,y,z). (2.4) 

i 

The current is given by 

J^^ = d^ = e ^(ii)^ = e J^ivi)^. (2.5) 

i i 
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The expectation value of the current is 



I 

i 

h 



— oo 



di'Tre-^^°(*-*')/'^[//'(t'),Po]e'^°(*-*')/^J^ 
We assume the time dependence of as {Juit)) — J^{u;)e~'^'^^, then we 

where the conductivity is written as 



/■oo 

Jo 



Here we have defined 



Due to the relation 



[a„, po] ^ -ihpo dXa„{-ihX) , 
Jo 



we obtain 



Jo 

= / dX{j,{-ihX)j^{t)). 

Jo 



Since the time-derivative of d> is written as 



4- = ~Mt)J4o) - J.(o)j^(t)), 



the conductivity is given by 

/■oo /-c 
^0 -'0 



oo ^i(uj+i8)t 

dt- 



i(a; + i5) 

/oo ^i{to+iS)t _ 



where 



QUt)^-fm[Mt),Mo)]), (2.14) 

and Q^yijjj) is its Fourier transform. Q^^{uj) is evaluated from the analytic continuation of 
the thermal Green's function: 

Q?.!^) = Q^.S^n ^nu + iS), (2.15) 

Q^^ir) = -(TJ^(r) J,(0)) = ^ J] gM.(^u;n)e-'""^ (2.16) 

Prom these equations we can derive the sum rule for a/j^i,. Let us define the Fourier transform 
of 0^i.(i) as 

/oo 
dt<l>^,{t)e'^\ (2.17) 
-oo 

then an^{u!) is written as 

Hence the following formulae are followed: 

/oo 
di^a^^icu) = TTcf>^,{t = 0), (2.19) 
■oo 



\imuja^,{iu) = i(l)^,{t = 0). (2.20) 
Since = (l/m)(pj^ - (e/c)A^) and = eX)i('yj)M> we obtain 

(Pnu{t = 0) = -Tr[a^,po]^M 

3 

= (2.21) 



Then the sum rule is written as 



/ (ia;Rea(a;) ^"Ltll^ (2.22) 
Jo 2 m 

In this derivation the translational invariance of the potential term is important since we 
used the relation xi^ — Vi^ = (l/m)(pi^ — {e/c)A^). In the Drude formula 

a{u;) = (2.23) 
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the sum rule is clearly satisfied. 

Please note that the above formulas are derived for the uniform external fields. An 
extension to the spatially oscillating fields is performed in a straightforward way. Here we 
set ^ = 1. Let us consider the spatially varying applied fields: 

H'{t)^^- [ drUr)E,{r,t), (2.24) 



c 



or 



where we assume 



H'it)^-^ J drj,{v)A^{r,t), (2.25) 



E^{r,t) = e^e'^^-'-^'\ (2.26) 

and A^(r, t) is the vector potential satisfying En{r, t) — (— l/c)^^(r, t) — {iu;/c)A^{r, t) and 
divA = 0. We set ^^^(r, t) — e J2i ^i/i<^(i" ~ ^i) ^-^d J/* is the current operator given by 

i 

The conductivity is defined as the coefficient of the hnear response of the current to apphed 
fields: 

J^(r,i) = (7^^(q,a;)£;^(r,i) 

= fdr' f dt'a^,{r -v';t- t')E,{r', t'). (2.28) 

J J —oo 

The expectation value of j//(r) to the first order in H' is evaluated as 

= i P dt' j drIVe-^^°*'[e,(r),po]e'''°*'e.e^'>"-e-^'^(*-*')j^(r) 

= -y ciiV'^*'Tre-*^°*'[e.(q),Po]e*''°*'e.e-^'^%(r), (2-29) 

where ^^(q) is the Fourier transform of ^^^(r). We follow the same procedure as for the 
uniform external fields and note that ^//(r) = J^(r), then we have 

<P^M:t)= dX{Mq,-thX)j^{r,t))e-''^-''. (2.30) 
Jo 

Here we neglect the A-term in the current since this is the higher-order term in external 
fields and we take a spatial average to write 

(q) = / drj,{r)e^^-^ = ^ ^(p.e^'*-^' + e^'^-'-'p,)^, (2.31) 

i 
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0^,(q,i)=/ dX{Uci,-ih\)j^{q,t)), (2.32) 

^0 

POO 

(J^M^u)^ dte'^'(t>^M,t). (2.33) 

^0 

Now let us define the retarded response function: 

X^.(q,i-0 = -i^(i-0([jKq,^),>(q,i')]), (2-34) 
and its Fourier transform given by 

/oo 
rfte*'^*i^^,(q,i). (2.35) 
-co 

The formula for the optical conductivity is followed as 

(7/.^(q,^) = ^^^^ [K^u{,(i,uj) - K^^{(i,Q)]. (2.36) 

The current response function is written in the form, 

( , ■s:\ v^e"^-^" (n|j,,(q)|m)(m|jt(q)|n) (n|j;(q)|m)(m|j^(q)|n) 

nm 

where |n) denotes a complete set of exact eigenstates of with eigenvalues and Z is 
the partition function Z = X^„e~'^^". The imaginary part of i^^j/ is given by 

ImX^,(q,a; + i5) = -7r(l - e"^-) J] ^^5(^„ - - u;)(n|jt(q)|m)(m|j,(q)|n). (2.38) 

nm 

The retarded Green's function K^i,{q,uj) is evaluated from the thermal Green's function 
through the analytic continuation: 

K,^u{c[,i^) = Kf,^{q,iUn ^ uj + i5), (2.39) 

where 

[0 . 

K^,{q,iujn)= dTe'^-^K^4q,T), (2.40) 

^0 

K,M:r) = -(riJ(q,T)j.(q,0)). (2.41) 
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In order to calculate the conductivity for the uniform external fields, we take the limit 
q — > 0. For the direct-current conductivity, we must take the limit q ^ first before uj —>■ 0. 
The current operator is given by 

j^^) = ^ Z1(P + ^q)4+q-Cpa = ^Y1 P4+q/2,aCp-q/2,a. (2.42) 
per pa 



In the limit q /^^(^(q, = 0) is evaluated as 
2e^l 



n k 



n k n k 

5^^, (2.43) 



a^,(cj) = ^— [/r^,(q= 0,cu + i5) + 6^,]. (2.44) 



m 

where Go(k, ie„) is the Green's function for the non-interacting electrons: (jo(k, ie„) ■ 

{ien — ek)~^ for e„ = (2r;, + l)7r//5. Then the uniform conductivity is given by the formula 

i Ne^ 
— -[K ^,{ci=0, 10 + i5) + — 
u + id m 

The current response function K^y for q 7^ in the non-intracting system is 

( ■ \ 26^-^ /(ep-q/2) - /(ep+q/2) ,^ 
Ki,y{q,tuJn) = 5- yp^Pu— ■ . (2.45) 

p + - ep+q/2 

Since -ft'^;^(q, u) vanishes in the limit q ^ for finite 7^ 0, we have the sum rule 

dujReaf,u{uj) = 6/,^, (2.46) 

Q Zi in 

and the Drude weight 

D = (2.47) 

m 

as the coefficient of the delta function: Yiea^yiijj) = D5{ijj)5^u- The sum rule for q 7^ is 
derived similarly. The commutator as in eq. (j2.21|) leads to 

I du;Recr/,,,(q,u;) = 7rRe0/,,,(q,t = 0) = 6^y. (2.48) 

J -00 ^ 

B. Mattis-Bardeen Formula 

The infrared absorption formula in BCS model was first derived by Mattis-Bardeen. 
We use the standard notations for the Green's functions; 

G(p,r) = -(cp.(r)ct (0)), (2.49) 



F(p,r) = (c_pi(r)cpt(0)), 



(2.50) 



F\p,r) = {cl,{T)cl^,m, (2.51) 
and their Fourier transforms given by 

G(p, ze„) = + (2.52) 

F(p, ie„) = F^(p, ie„) = -u^v^{- - - — ^-^), (2.53) 



where E^ = + , = (1/2)(1 +^p/_E'p) and Vp = l — u^ for ^p = Cp — /i. We evaluate 
the current response function written as 

2g2 ]^ 

Kf,^{q,iu£) = —y2pf,Pu-^y2[G{p-ci/2,ien)G{p + q/2,ien-iuJe) 

p n 

+ F{p - q/2, te^)F\p + q/2, ze„ + iue)]. (2.54) 
We set ia;£ —iuj£ in the first term, then using the symmetry q ^ — q we have 



lit j-i \ 



Ep^Ep_ Ep^Ep_ J iiOi + Ep^ — Ep 



+ ifiEp,) + f{Ep_)-l)- + 



4 \ Ep^ Ep_ Ep^Ep^ J iiOi + Ep^ + Ep_ 

+ (1 - /(i^pj - /(i^pj)7{(i - + #-) - #4^}- 



4 -Ep^. -Ep. -Ep^. -Ep^ iuJe — Ep^ — Ep_ 

(2.55) 

where p+ = p + q/2 and p_ = p — q/2. Here we consider (i) dirty superconductors or 
(ii) thin films satisfying d <^ ^ ^ vp/A for the film thickness d and the coherence length 
^. In these cases we can regard p+ and p_ as independent variables. In the case (ii), since 
qvp ^ A ~ we can do Abrikosov's replacement |l| , 

5^-iV(0)^ / dUd^p_. (2.56) 

Then in the isotropic case we obtain the Mattis-Bardeen formulae for ai{uj) = Iiea{uj) = 
— {l/uj)lmK{LJ + i6) and (72{uj) = Imcr(co') = (l/ci;)Re-ft'(ti; + i6): 



1 puj-A / a2 

^{co) = - / dEN{E)Nicu -E){1- 2f{u + E)) (l - — ) Q{uj - 2A) 

fin J A \ t,[UJ - t,) ^ 





/•cj— A 


w J 


A 


1 




+ 2- 






Jo 



A' 



E{uj + E) 



(2.57) 



where (Ti„ is the real part of the conductivity for normal state. In the above formula for (72, 
the integral is calculated as follows. 



4 EpEpi u — Ep — Ept 

= / de dEdE'N(E)N(E')(l )P 6{e + E - E') 

J-oo J A EE' u-e 

l-OO l-OO A 2 1 

= / de dEN{E)N{E + e){l-—^—-)F 

Ja-e J a E{E + e) u-e 

/oo poo 
de / dEN{E)N{E + e){l- , 

/■oo ^2 

^TTlmy^ dEN{E)N{E-u^){l-^^^—^) 

J rrii 



-1.1 1 



E(a; - - A^ 



' mox(aj— A,A) 

= TT / dF, ^ ' 2.59 

imax(A-c.,-A) ^/ A"" - E'' ^ {uj - Ef - /S.'' 

At T = 0, the integral has the form for complete elliptic integrals by changing variables of 
integration to x where x — {2E — uj)/{ijj — 2A) for (Ji^; 



1 ^ Ja 



(Tin C^Ja - AV(^ - Ef - A2 

/ , 2.60 
io y(l-a:2)(l- 



where k = {uj — 2A)/(a; + 2A) and u) > 2A. We then obtain the Mattis-Bardeen formula: 

— H = ((l + — ^ E{k) - —K{k)] e{u - 2A), (2.61) 
(7i„ \\ o; y a; / 

where 



m= I xlK^dx, K{k)= I -^=======dx, (2.62) 



Jo - x^){l - k^x^) 

are complete elliptic integrals. An expression for valid for all cu is 



£^H.(^4)^(,0.(§-i)i.(*', (2.63) 
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where k' = {1 - k'^f'^. 

In evaluations of K^y, the matrix notations are also employed for Green's functions in 
superconductors, 

G'(k,ze„) = (2.64) 

+ + Ak 

where Tj (i = 0, 1, 3) are Pauli matrices, and Ak are generalized to include the self-energy 
in the form: = Cn — So(e„) and Ak = Ak + Si(e„). The response function K^^^ is written 
as neglecting pair vertex corrections 

K^y{(l,iujm) = ^ V'p^.p^- y^TrG'(p+,ze„ + zcj„)G'(p_,ie„). (2.65) 

p n 

One can reproduce the Mattis-Bardeen formulae from this expression in a straightforward 
way. 



C. Optical Conductivity in London Superconductor 



In this section, we investigate the superconductor in the London limit: ^ ^ A for the 
coherence length ^ and the penetration depth A. Sin ce q£ -C 1 holds, we must consider the 
limit g — > for the response function in eq. 

K (n 1 - ""^^^^^S^h k ^ V + ''^'"^ + ^'^-^'^+ Ak_ Ak^ 



Sin ce q£ -C 1 holds, w 



We use the notation z = ie 



e^kl 



ien — iT,Q{ien) = z — zSo(z). Then we have 

z{z + iujm) + ^k + AkAk^ 



(2.66) 



(2.67) 

where A^^ = A + Si(z + icu) and C is the contour surrounding the poles of f{z) in the 
clockwise direction. We set z^ = z + iuJm, and write the integrand in the form, 

z~z+ + ^l + AlAl^ _ 1 



Z2 



;;2 



+ 



ZZ+ + A^A^^ + ~z: 



(e^ + A2-52)(e2 + A^ -(z + z^. 



(2.68) 
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Then the momentum summations are performed in the following way: 



/ X , 1 /" , r d zfiz) f(z) 
= NiO){- / dz[-- Jl^ + ii^ ])^, (2.69) 

where {■ ■ denotes the average over the Fermi surface. In the last equality, the first term 
gives only a constant contribution. The second term in eq. ()2.68|) gives 

1 1 



Z + lUn 



P~AiJ~zl~Ai^ 



X 



5i - AL 



(2.70) 



Then the current response function is written as 



2 1,2 

- ^N{0){k^K I dzfiz) 



X 



V 



c 



ZZ+ + AkAk 4 



2--Ai^~zl-Alj 



z^-Al + J~zl-Al 



(2.71) 



Let us consider the limit of weak impurity scattering, and write the renormalized fre- 
quency z and the superconducting gap Ak in the following forms, respectively: 



z = z + iTii 



Ak = Ak + zr2( 



P-Ai 



Ak^ 

P - Al 



-)k'- 



(2.72) 



(2.73) 



For the isotropic superconducting gap, we set -u = z/A and v = uA; then we have 

u 



v(z) = z + iT- 



z + iT- 



(2.74) 



V^F^ Vv^ - A2 

where F = Fi — F2. The relation \/ P — A"^ = ^/v'^ — A"^ + iT2 is also followed. The density 
of states is given by 



Ns{z) = N{0)Re 



u 



(2.75) 
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In the limit as F — 0, this reduces to 

N,{u) = N{0)^^L= for 1^1 > A, (2.76) 
and Ns{uj) = otherwise. For the anisotropic case satisfying (Ak)^ = 0, v{z) satisfies 

v{z) = ~z = z + ^T{^^^===)y (2.77) 

for F = Fi (F2 = 0). Let us examine this case in more detail; K^^ is given by (where we use 
the notation v+ = v{z + iuJm)), 



X II ^^^ + ^1 V, 

^-iV(0)(A;^A;, / dzfiz' 



Now V has a cut along the real axis; as the cut is crossed, the continuation is performed as 
follows 0, 

V — > V*, 

-(^^=Y. (2.79) 

The complex integration is reduced to integrations along Imz = and Imz = uOm- We obtain 
the expression for the current response function from the analytic continuation: 



„2l2 rco 

K(q = Q,u + i5) = ^NiO) J de[{f{e) -fie + co)) 



r 1 vje + u) vje)* 

\ / 7 '. To I 7 To *~o"/ fcJ 



v{t^u)-v{ty^v{e^uf-Al V^(e)2- A2'' 
- /(e + c^) 

X r 1 t;(6 + a;)* t;(e)* 

^^;(e + a;)* - ^;(e)* ^^(e + cu)*2 - ^t;(e)*2 - A^ 

, 1 / ^(e + ^) v{e) 

^ z;(6 + c.)-t;(e)^y^;(7T^J)^3A^ ' 
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where the average K^x + Kyy is written as 2K (in two dimensions) and we use v{uj — iS) 
v{u! + i6)* . If we use the relation in eq. ()2.77p . the expression for Imi^' is simphfied as 

lmK{q = 0,u; + zS) = ^iV(0)|, rfe (^tanh (^^^ - tanh + 



X Re ^ - ^ -— . (2.81) 

\v{e + uj) - v{e) v{e + uj) - v{e)* J ^ ' 

In the collision less limit F ^ 0, an expansion in terms of F gives the conductivity, 



X (Re ' ' =)r, (2.82) 

where 0\n = {ne^r) /m/ {ut)"^ ^ {ne'^T)/m ■ l/[(a;r)^ + 1] with r = 1/(2F), and we use the 
density of states in eq. ()2.76p in the limit F — > 0. 

For the rf-wave symmetric superconducting gap in two dimensions, the average over the 
Fermi surface is given by the Elliptic function, 

1' I r'^'" V 2 A 

--)k = 7^ # . . ... =^ = (2.83) 



A2cos2(20) ^ 27r7o Vt;^ _ ^2^052(20) 



TC V 



for V > and = Acos(20). Since the relation K{l/x) = —ixjyX — x^- K{\j \J\ — x^) 
holds for X > 1, the renormalized frequency u[lS) = v{uj)/A is written as 

u .T , u 

= I + r 2 „ ^ 1 

We will do a continuation of the elliptic integral to a general complex argument, the real 
part of the optical conductivity for the (i-wave superconductor in the London limit is given 
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by[ll| 

2^ 



)— / cte- tann — — tann 

ujLoo 2V V2; 



X im(-z)-{-— —-——-[— =^===K 



Mfe)* / 1 



*2 



1 2 , n(e + cu) / 1 



(2.85) 



The imaginary part Ois is also obtained from eq. ()2.80p as 

G2s{^) = Ima,(cj) = -ReK{q = 0,uj + i6), (2.86) 
where Rei^r(g = 0, a; + i6) is written as 

ReK{q = 0,u; + z6) = ^N{0) j de\U{e) - /(e + u)) 

X Re(- J -— , ^ 

1 



v{e + uo) — v{e) t: v{e + uo) v{t) 
Vie + a;)Re— — i ^ "(^^T^) " ^^^Y^)] 



^Ar(O) j ^rfe[-(tanh(^) -tanh(-^^ ^ 



)) 



, .1 1 2 u{e + u) ^. 1 



- , "^^^^ K{ ' )) 
1 1 2, u(e + cj) 1 



X 1 1 2, uie + uj) , 1 

- 2/ e + cu Im-- ^ , K{— 

A u{e + iu)- u{e) tt y 1 - „(e + uj)^ ' ^1 - u{e + ujy 

""^"^ -M^=^—))). (2.87) 
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D. Conductivity Sum Rule 



As shown in Section II. A, the sum rule holds for the conductivity: 

2 

71 TIC 

Rea{iu)duj = , (2.88) 

2 m 

where cr{uj) is divided by the volume so that the quantity is of the order of 0(1) and n = N/V 
is the electron density. In superconductors, there is a dramatic change in the optical con- 
ductivity stemmed from opening of an excitation gap. The change of the conductivity is 
compensated by the formation of a zero frequency 6 function peak to preserve the con- 
ductivity sum rule.jssi From the general formula in eq. ()2.36p . the real part of the optical 
conductivity is written as 

Recr(a;) = -lm—^—^{K {q,uj + i6) - K{q,uj ^ 0)) 

= --lm(K(q, u; + i6)- K(q, u ^ 0)) + n6(uj)Re(K(q, u + i6) - K(q, u 0)) 

OJ 

= - -lm( K (q,uj + i6) -K(q,uj ^0)) + 7i6(uj)ojlma(u). (2.89) 

OJ 

In superconductors, the first term indicates ais = Rea in the finite frequency region whose 
weight is removed by the opening of the gap, and the second term represents the condensate 
peak which recovers the lost weight for finite oj. The Drude weight is then given hj D = 
nlimi^^oojlmcr^oj). The sum rule is expressed as 

Rea{oj)doj = - -lm{K{q,oj + i5) - K{q,oj ^ O))doj + - \imojlma{oj). (2.90) 

Jq OJ 2 t^^o 

If the total weight is conserved in the superconducting transition, we have the relation from 

D TT 

{ain{uj) - ais{uj))duj = — = - lim uja2s{uj), (2.91) 

where denotes the imaginary part of a in the superconducting state. In superconductors, 
the weight of the condensate peak can be estimated from the measurements of the London 
penetration depth, 

\imoja2s{(^) = —To = (2.92) 

where the right-hand side is written as l/(A|/io) in MKSA unit using the permeability 
of vacuum /iQ. As we will show in the next Section, the sum rule in eq. ()2.91|) actually 
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holds for the optical conductivity spectra cr{uj) of NhNi-^-C^^ obtained using the Refiectance- 
Transmittance method. 0] 

Recently, however, in cuprate superconductors, there is an experimental report that the 
sum rule is violated for c-axis conductivity. [■3] The weight estimated from the reflectance data 
is much larger than can be accounted by the spectra integration in eq. (j2.91|) . This issue is 
not resolved at present whether the extra weight is coming from outside the frequency range 



4l| 



measured or coming from the lack of quasiparticle poles in the Green's functions. 

At the last of this Section, we briefly discuss the ration between the conductivity sum 
rule and the f-sum rule. The dielectric function e(u;) is related to the conductivity through 



the relation 



42 



e{uj) = l-—a{uj). (2.93) 



The f-sum rule states that 



eH ^ 1 - ^ (iu^oo), (2.94) 



UJ 



where Up = A'nne^ j m is the plasma frequency. This indicates 



o{uj) (w ^ oo), (2.95) 

ijj m 

and hence in the high-frequency limit lma{uj) = {ne^ / m) / uj . Since e{uj) is analytic in the 
upper half-plane, performing an integration of uj{t{uj) — 1)) in the upper-half plane, we obtain 

dujiA'Ka{uj) = incUp. (2.96) 
This equation implies the sum rule in eq. ()2.88|) . 

III. REFLECTANCE-TRANSMITTANCE METHOD 
A. Method of Analysis 

The conventional FIR spectroscopy based on a Kramers-Kronig (K-K) transformation, 
however, is rather unfavorable for studying electronic properties in small energy region. 



52| | Recently, a new method to examine far-infrared (FIR) spectroscopy has been developed 



without devoting to evaluating the Kramers-Kronig transformation. 27| In this method the 



optical conductivity is estimated from the data of reflectance spectra R{ll!) and transmittance 
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spectra T{u) by substituting them into a set of coupled equations. The new method is free 
from the conventional difficulties in far-infrared reerion since we do not need the aid of 
Kramers-Kronig transformation. This method is referred to as the R-T method since 
both R{uj) and T{uj) are necessary for analyses. The basic concept of the R-T method was 
first applied to the study of <j{uj) of superconducting NbN thin films deposited on MgO 



and Si substrate for u below ~ 150cm~^ at T = 5 — 20K. 26|, In this section we briefly 
describe the concept of the R-T method and its advantages over the conventional method 
based on K-K analysis using bulk samples. 

The reflectance spectrum Ri{uj) and transmittance spectrum Ti{u) of a single- layered 
system (such as MgO substrate) are given by 

R,{uj) = \n{uj)\^, (3.1) 



T^iio) = \h{uj)\\ (3.2) 

when the incident radiation is introduced in the direction normal to the layer surface. Here 
ri{uj) and ti{uj) are complex reflectance and transmittance, respectively. If the multiple 
internal reflections within the layer are exactly taken into account, ri{uj) and ti{uj) are given 
by 

(1 - A^i)(iVi + l)e^('^/'=)^i'^i -f (1 + iVi)(iVi - l)e-^('^/^)^i'^i 
""^^"^^ ~ ~ (1 - Ni){Ni - l)ei(-/c)N,d, + (1 + jv^)(jv^ + i)e-i(a./c)iVidi ' (^.3) 

4jY^g-i(<^/c)di 

^'^"^^ ~ ~ (1 - Ni){Ni - l)e^(-A)^i'^i + (1 + Ni){Ni + l)e-^(-A)^i'ii ' ^ 
where c is the velocity of light in vacuum; di is the thickness of the layer; Ni = ni + iki is 
the complex refractive index. For a two-layered system (such as YBCO thin fllms deposited 
on MgO substrate) placed in vacuum, the reflectance R2{uj) and transmittance T2{u!) are 
given by 

R^icu) = \r2{uj)\\ (3.5) 

T2(cu) = |^2(a;)|^ (3.6) 

where we assume the same conditions for incident radiation. r2{uj) and t2{u!) are given by 

A+B+C+D 
^^^^^ = E + F + G + H^ ^'-^^ 
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where 

B = -{Ni + N2){N2 - l)(iVi + l)e-^(Wc)(Af2d2+iVidi)^ ^g^^g^ 

C={Ni + N2){N2 + l){Ni - iyi-^/c)(N2d2+N^d^) ^ ^3_^^^) 

D = {Ni- N2)iN2 - l)iNi - i)Q-ii-^/c)(N2d2-N,d,) ^ ^3_^2) 

E = {Ni- N2){N2 - l)(iVi + i)ei(Wc)(Af2rf2-Aridi)^ (3^^3^ 

F = {Ni + N2){N2 + l)(iVi + i)e-i(-'/c){N2d2+mdr)^ ^3_^4^ 

G = -{Ni + N2){N2 - l){Ni - i)e^(^/'=)(^2rf2+7Virfi)^ (3J5) 

H = -{Ni - N2){N2 + l){Ni - i)e-*(^/^)(^2rf2-Afidi)^ (3^^g) 

J = 8NiN2e-'^'^/'''^^'^^+'^'\ (3.17) 
Here d2 and iVo = 712 + ik2 are the thickness and the complex refractive index of thin films, 



respectively. j2^, |22| The equations ()3.1|1 - ()3.17|1 are basic relations in R-T method. 

The coupled equations are numerically solved to determine the values of complex refrac- 
tive indices rii and fcj {i = 1,2) as functions of uj. First, the coupled equations for Ti{uj) 
and Ri{uj) measured for the MgO substrate are solved using the Newton method, then we 
obtain rii and ki as a function of uj. Second, T{uj) and R{uj) of YBCO/MgO are measured. 
The measured values are substituted into R2 and T2 for a given value, as well as rii and 
ki for the MgO substrate obtained in the first procedure. Then n2 and k2 of YBCO are 
determined as a function of uj by solving the coupled equations numerically. 
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Here, we briefly discuss both the advantages and disadvantages of the R-T method in 
comparison with the conventional method based on the K-K analysis using bulk samples. 
First, we can estimate precisely the optical constants of materials such as the complex 
refractive index. Second, the accuracy of cxi obtained by the R-T method is expected to be 
better than that obtained using the conventional method in the low frequency region. We 
have three reasons for this expectation, (i) In the conventional method we use only R{uj), 
while in the R-T method both R{uj) and T{uj) are employed to obtain the complex refractive 
index, (ii) We must be careful about the K-K analysis in the conventional method. For 
metallic materials, R{uj) takes high values in the small u region and thus the experimental 
errors increase through the K-K transformation. The errors in R{uj) propagate to the errors 
in (Ti. The conventional method is not so reliable in the study of (Ji[uj) in the small u 
region for metallic materials, (iii) The K-K analysis usually requires an extrapolation of 
R{uj) to the low- and high-frequency regions beyond the measurable u region. Since there 
are no guiding principles for the extrapolation of R{uj), there may be some uncertainty in 
the K-K transformation in the conventional method. In particular, the ambiguities in the 
extrapolation are not negligible in the low frequency region. The R-T method does not have 
this kind of difficulty. 

A disadvantage of the R-T method lies in the fact that the range of u for which the method 
is applicable is limited because the substrate must be transparent to the incident radiation 
to measure T{uj). For MgO as the substrate, the transparent range is about — 300 cm~^ in 



far- infrared region even at T below 10 K. [2^123 In addition, the high-frequency limit of the 
transparent range decreases steeply as T increases for MgO; this limit is approximately 100 
cm~^ at room temperature. Since the conventional method based on the K-K transformation 
is expected to work very well for u > 300cm~^, the region of uj within the scope of the R-T 
method is < 300cm~^. 



B. Experimental Results 

In order to examine the feasibility of the R-T method, we report R{uj) and T{uj) of 
NbNi_a;Cx thin films deposited on MgO substrates in the far-infrared region. The substrates 
were 0.5 mm thick. The thickness of the NbNi.^C^. layers was about 40 nm. The films were 
epitaxial and Tc was estimated to be 17.5 K by electrical measurements. The value of x was 
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FIG. 1: Reflectance spectra R{u>) of NhNi-^^C^: tliin films deposited on MgO substrates at T = 4.3 
K (solid line) and T = 20 K (dotted line). 




FIG. 2: Transmittance spectra T{uj) of NbNi_a;Cj; thin films deposited on MgO substrates at 
T = 4.3 K (solid line) and 20 K (dotted line). 
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FIG. 3: (Ti(cj) for NbNi_j;Ca; calculated by the R-T method at T = 4.3 K (solid circles) and 
T = 20 K (open circles) . The dashed line shows the results of calculations using the Drude formula 
at T = 20 K. 
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FIG. 4: a2{io) for NbNi_^Ca; calculated by the R-T method at T = 4.3 K (solid circles) and 
T = 20 K (open circles) . The dashed line shows the results of calculations using the Drude formula 
at T = 20 K. 
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FIG. 5: The relative conductivity ratio ais{io) / ain{uj) where cri(uj) at T = 4.3 K is divided by 
(Ti(u;) at T = 20 K. The sohd hne shows the results obtaned from the Mattis-Bardeen formula. 
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FIG. 6: The relative conductivity ratio cT2s(w)/cr2n(w) where a2{uj) at T = 4.3 K is divided by 
(72 (w) at r = 20 K. The solid line shows the results obtaned from the Mattis-Bardeen formula. 
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estimated to be less than 0.3. 

The electrical resistivity measured by the four-probe method was p = 5.2 x 10~^ Qcm at 
T = 20 K. The carrier density in the normal state n was measured by the van der Pauw 
method: n = 1.29 x 10^^ cm~'^ which is slightly lower than the value of = 2.39 x 10^'^ cm~^^ 
reported previously for NbNjs^]. The carrier scattering rate F at T = 20 K was estimated 
as r = pe^n/m* = 1.9 x 10^^ s~^ = 6.3x10^ cm~^, where the effective mass m* is assumed to 
be equal to the electron rest mass. The superconducting gap has been estimated as 2 A ^ 50 
cm^^. Thus NbNi_a,Ca; is suggested to be a typical dirty-limit BCS superconductor. 

R{uj) and T{uj) obtained at T = 4.3 K and 20 K are shown in Figs. Q and El The 
interference fringes due to multiple internal reflections within the MgO substrate are clearly 
visible in both figures because the MgO substrate is highly transparent in this u region at 
T = 4.3 K and 30 K, and because the NbNi-i^Cj; film is thin enough to transmit far-infrared 
radiation. R{uj) at T = 4.3 K exhibits an obvious reflectance edge at ~ 65 cm~^ and 
R{uj) ~ 1 for u; less than the reflectance edge frequency, which is a special characteristic for 
superconductors. T{uj) at T = 4.3 K exhibits a maximum at ~ 60 cm^^, which is related 
to the evolution of the reflectance edge in R^u). 

We show cTi^u) and cr2{uj) spectra Figs. El and HI for NbNi_a.C2, calculated by the R- 
T method using the experimental results shown in Figs. ^ and The value of the dc 
conductivity cri(O) at T = 20 K is estimated to be cri(O) ~ 2.0 x 10^ i7^^cm~^ from Fig. 
13 this value agrees well with the value of 1/p = 1.9 x 10^ Q~^cm~^ estimated from the 
electrical measurements. 

We evaluated the relative conductivity ratio ctis^lj) / ain{uj) and a2s{uj)/ain{uj) from the 
results in Figs. El and (H where ain is ai at T = 20 K and ais and are at T = 4.3 
K, respectively. The real and imaginary parts of the relative conductivity ratios are shown 
in Figs. El and respectively. Here theoretical curves obtained using the Mattis-Bardeen 
theory are also shown by solid lines, where we set 2A = 52 cm~^ in accord with the value 
reported by the junction method (28l|. The ratio of 2A to Tc is given by 2A/ ksTc ~ 4.3, 
suggesting the strong-coupling superconductivity in NbNi_2.Ca;. 

The experimental results for (Jis/ain in Fig. El exhibit an excellent agreement with the 
Mattis-Bardeen theory, and cr2s/o"i„ also shows a good agreement for u less than ~ 60 cm~^ 
as shown in Fig. El The agreement is, however, poor for uj larger than ~ 70 cm~^ in Fig. EJ 
this anomalous behavior may be due to impurities 66 1. 
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Now we investigate the conductivity sum rule in eq. ()2.91|) . From the integration of spectra 
(Tin and ais, Xl was estimated to be ~ 193 nm, while the a2s spectra in the superconducting 
state gives Xl ^ 200 nm. These values show an excellent agreement, and also agrees well 



with the value reported previously. 



28l | This indicates that the sum rule holds for NhNi^^Ca,. 



IV. ELECTRON-DOPED HIGH-T^ SUPERCONDUCTOR: LONDON LIMIT 



Oxide high-Tc superconductors have been investigated intensively over the last decade. 
The (i-wave superconductivity is well established for hole-doped superconductors. However, 
there is a class of high-T^ superconductors doped with electrons, for which both 

s-wavej3| and rf-wave pairing^, 13, have been reported. Nd2-a;Ce2^Cu04 is a typical 
example of electron-doped materials and the symmetry of Cooper pairs has been contro- 
versial. It is important to examine the symmetry of Cooper pairs in the study of high-Tc 
superconductors . 

Since the superconducting gap A in Nd2--xCe^.Cu04 is very small, there have been no 
reports on the study of the nature of the superconducting gap of Nd2_xCe2:Cu04 through 
such techniques, although there have been a number of reports on the FIR spectroscopy of 
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FIG. 7: Real part of the optical conductivity as a function of oj for several values of temperature. 
From the top T/A = 1/(3 = 1, 1/2, 1/4 and 0. 
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FIG. 8: Optical conductivity calculated from the R-T method at T = 4.3K (circles) and 30K 
(squares). The diamond indicates the DC value at T = 30 K. 
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FIG. 9: Optical conductivity (circles) by the R-T method and theoretical predictions at T = 
(solid curves). Form the left 2A = 50cm~^, 60cm~^ and 70cm~^ 
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Nd2_a.Ce:,.Cu04.|2ll|43| 

The purpose of this paper is to investigate FIR optical properties of Nd2-xCea;Cu04 
obtained by the R-T method from a viewpoint of unconventional superconductors. We 
will show that the available data for the optical conductivity and transmittance are well 
explained by rf-wave pairing model in the clean limit. The value of superconducting gap is 
estimated as 2 A ~ 60 — 70 cm~^, which is consistent with the available value estimated by 
.ca„.,,tu„.e.,„gspe*.o.oop,y 

The frequency dependent conductivity cr{uj) was calculated by Mattis and BardeeUjjs^ 
Abrikosov et al. \2\ and Skalski et al. j^l for isotropic superconductors. The original Mattis- 
Bardeen theory was carried through for a conventional type-I s-wave superconductor, where 
the coherence length ^ and magnetic penetration depth A satisfy ^ ^ A. The opposite 
limit ^ <C A (London limit) was also examined for s-wave pairing by field theoretical 



treatments. 



45| For the high-T^ compounds of type-II superconductor with small coherence 



length, the formula in the London limit is appropriate for optical conductivity measurements. 
Recently the conductivity <y{uj) of an unconventional superconductor has been examined the- 



3.0 




FIG. 10: Observed Transmittance and the theoretical curves at T = (solid curve) for 2A = 
60cm~^. 
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oretically in the London limit, 
in Section II: 
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mm 

We use the current response function shown 



K^iy{q, iujm) = — f 22 ^fj-^i^lS z2 '^^[^('^+' + ^^m)G'(k_, ze„)], (4.1) 



k ^ 



where k± = k ± q/2 and = {2n + l)7r//5. The single-particle matrix Green's function is 

^ ' (^n-S(6„))^ + ei + A^ ' ^ ^ ^ 

where is the anisotropic order parameter and S(e„) is the self-energy due to impurity 
scattering, r* (z = 0, 1, ■ ■ ■) denote Pauli matrices. Since we consider the case where ^ <C A 
holds, the real part of optical conductivity is well approximated by the formula in the London 
limit: 

cris,f,u{uj) = --\imlmK^^{q,Lj). (4.3) 

Our focus is the collision less limit of the normalized conductivity to compare it with the 
data for Nd2_i.Cea;Cu04 since ^ -C £ holds for the mean-free path i. For anisotropic su- 
perconducting order parameter such that the average over the Fermi surface vanishes 
(Afc) = 0, the expression for au = cris,xx{^^) in the collision less limit on the plane is simply 
given by 19| 

-7^ rfx(Re J-±-= ){Re ^ ^ ==)[tanh(^) - tanh( ^^ ^ 



(4.4) 

which is an angle-dependent generalization of the Mattis-Bardeen formula. For the d-wave 
symmetry, the average over the Fermi surface denoted by the angular brackets is defined as 

(Re , ) = Re / ^ , (4.5) 

^/x^-Al' J 2tt ^x2-(Acos(20))2' 

where the order parameter is factorized as A^ = Acos(20). In FigQwe show the behaviors 
of (Tis/crin as a function of u for several values of temperature T. The infrared behaviors 
reflect the lines of nodes on the Fermi surface. 

FIR reflection R{uj) and transmission T{uj) measurements were performed for 
Nd2_a.Ce2,.Cu04 {x = 0.15) thin films deposited by laser ablation onto (001) MgO sub- 
strates. The thickness of Nd2-a;Cea;Cu04 thin film was about 40 nm. Tc was estimated to be 
~20K. The electric field of the FIR radiation was predominantly parallel to the a-b plane. 
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The conductivity spectra were evaluated by the R-T method from the data for R{uj) and 
T{u) at T = 4.3 and 30K. 4^ 

The R-T method provides us rehable data of spectroscopy in far-infrared region for which 
comparison between the experimental data and theoretical analysis is possible. In the R-T 
method both the reflectance spectra R{uj) and the transmittance spectra T{uj) are measured 
experimentally from which a set of coupled equations are followed describing the transmit- 
tance and reflectance of a thin film on a substrate. The coupled equations are solved nu- 
merically by the Newton method to determine the optical conductivity. This method is 
free from the difficulties in the infrared region which occur commonly in the conventional 
method employing a Kramers-Kronig transformation. In FiglHl we show the real part of the 
optical conductivity obtained from the R-T method at T = 3.4K and T = 30K. In FigEl 
we show the observed data and theoretical curves at T = for 2 A = 60 and 70 cm~^. The 
experimental data a{3.4:K)/a{30K) normalized by the normal state values at T = 30K are 
shown in FigEl It is obvious from the experimental results that there is no evidence of a 




CO/A. 



FIG. 11: Transmission Tg/T/v for anisotropic s-wave models. The solid curve without marks is 
the Mattis-Bardeen result. Squares, triangles and circles are for the prolate (a = 0.5), ab-plane 
anisotropic (c = 0.5), and oblate (a = —0.5) gaps, respectively. The oblate form shows a small 
increase compared to other types. For the oblate gap, Aq = Amax, and for other types, Aq = A. 
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true gap, which is suggestive of anisotropic superconducting gap, since the spectral weight of 
conductivity should vanish for uj < 2A at T = in conventional isotropic superconductors. 
It is also shown in FiglHlthat they are well fitted by the curve with 2 A = 60 cm^^, which is 
consistent with the value estimated by scanning tunneling spectroscopy measurements, j^] 
Transmission curve is also presented in FiglTUl where Ts/Tj^, the ratio of the transmission 




FIG. 12: Transmission Ts/T^ for the two-band model. The data points (large solid circles) are 
taken from ref. 2^ at T = 6K. The solid line without marks shows the Mattis-Bardeen result with 
2A = 5meV. The dotted line is for the single gap of prolate type with 2/S.max = 9meV and a = 0.5. 
Triangles are for the singlc-gcip model of oblate form with '2,^Ynax 

= 9meV and a = —0.33. Others 

are for the two-band gap model where the ab-plane anisotropy for the c-band and the prolate form 
for the vr-band are assumed. The parameters are the following. Solid circles: 2/\raax = 8.5meV and 
c = 0.33 (fj-band: weight 0.45); 2 A 

max — 6meV and a — 0.33 (vr-band: weight 0.55). Open circles: 
2/S.rnax = 9meV and c = 0.33 (cr-band: weight 0.45); 2/S.max = 6meV and a = 0.33 (vr-band: weight 
0.55). Squares: 2/S.max = lOmeV and c = 0.5 (u-band: weight 0.4); 2IS.max = 7.5meV and a = 0.5 
(vr-band: weight 0.6). 
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in the superconducting to that in the normal state, is the experimentally measured quan- 
tity. The following expression for Ts/Tn is employed to determine the transmission curve 
theoretically, Q| 

= I (4 6) 

[TT + {^-TlJ'){ar/a^)]^ + [{l-T'J'){a,/ar.)f 

where ai and a2 are real and imaginary parts of the conductivity —{c/uj)K{ci, uj) for q — 0, 
respectively. Here we use the formula obtained from the two-fluid model for cr2. T/v is 
determined as Tjv — 0.05 from the expression for the ratio of the power transmitted with a 
film to that with no film given as 

= l/[l + and^]'. (4.7) 
n + 1 

Here d is the film thickness, n is the index of refraction of the substrate, and Zq is the 
impedance of free space. We have assigned the following values; d = 4 x 10~^cm, n = 3.13, 
Zq = 377Q and cr„ ~ lO^fi^^cm^^ and the Drude width is approximately equal to A cr„ is 
approximately given by the value at u = 0. Obviously the cj-dependence of measured trans- 
mittance agrees with the theoretical curve for 2A = 60cm^^ as shown in Figimi An agree- 
ment between the observed quantities and theoretical curve is remarkable, which should be 
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FIG. 13: Specific-heat-jump ratio to the BCS value. From the top the ratios for the oblate, 
o6-anisotropy, and prolate gaps are shown, respectively. 



31 



compared to the isotropic BCS prediction calculated from the Mattis-Bardeen equations. 







V. TWO-BAND ANISOTROPIC SUPERCONDUCTIVITY IN MAGNESIUM DI- 
BORIDE 



After the discovery of 39 K superconductivity in MgB2[39[, much attention has been 
focused on the study of its nature. An s-wave superconductivity (SC) was established by 



experiments such as coherence pea^ 



i in ^^B nuclear relaxation rate 



and its exponential 



dependence at low temperatures |35l. l65l|. An isotope effect has suggested phonon- mediated 

n 

s-wave superconductivity|8ii].In contrast to its standard properties, there have been several 
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FIG. 14: Real part of the optical conductivity for the two-band anisotropic model. The data points 
(large solid circles) are taken from ref.[2^ The parameters for solid circles are 2Amax = 8.5meV 
and c = 0.33 (a-band: weight 0.45); 2/S.max = 6meV and a = 0.33 (vr-band: weight 0.55). The 
parameters for open circles are 2lS.max = lOmeV and c = 0.5 (o"-band: weight 0.4); 2/S.max = 
7.5meV and a = 0.5 (vr-band: weight 0.6). The dashed line indicates the results obtained using 
the Mattis-Bardeen formula with 2A = 5meV. The dash-dotted line denotes the conductivity for 



the d-wave gap. 
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reports indicating unusual properties of the superconductivity of MgB2. Two different su- 
perconducting gaps have been reported: a gap much smaller than the expected BCS value 
and that comparable to the BCS value given by 2A = 3.53/crTI. Their ratio is estimated 



to be A„j„/A 

r, 



0.3 — 0.4 using several experiments 



49, 



551 . I65| . It is also 



reported that the specific-heat jump and the critical magnetic field are reduced compared to 



the s-wave BCS theory j6|, |65[. A strongly anisotropic upper critical field in c-axis-oriented 



1 



MgB2 films and single crystals of MgB2 is also observed j3|, 

The unusual properties of MgB2 suggest an anisotropic s-wave superconductivity or a two- 
band superconductivity. The band structure calculations predicted multibands originating 
from a{2pj.^y) and tt(2pz) bands. 30] In the ARPES measurements performed in single crystals 
of MgB2 three distinct dispersions approaching the Fermi energy were reported. jS^ | 



There have been several studies on the anisotropy of a superconducting gapiZ 

The two-gap model is shown to consistently describe the specific heat j7|, |40| and the 
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FIG. 15: Thermodynamic critical magnetic field Hc(T) normalized by Hc(T = 0)bcs- The bold 
dashed curve indicates data from ref. y| and the bold solid curve indicates those obtained using the 
present two-band anisotropic model. The thin solid curve indicates the BCS results. The results 
for the cr-band (ab-plane anisotropic) and the vr-band (prolate form) are also shown. 
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TABLE I: Anisotropic parameters in the SC gap function used to fit several physical quantities. 
The upper four rows are for the single-SC gap model and the last row is for the two-band anisotropic 
model for comparison. The cross indicates that we cannot fit experimental data by the correspond- 
ing z factor. A-i^ column //c2 indicates that experiments are explained qualitatively but not 
quantitatively, jjj] The big circle show that we can fit the data using the same parameters in the 
column The effect of a-band anisotropy is small. 
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upper critical field if^oj38l| with the adoption of the effective mass approach. 

Here we show that this material is described by two order parameters attached to a- and 
TT-bands. Two order parameters further must have different anisotropy to explain the ex- 
perimental results consistently. In this paper, we examine optical properties and thermody- 
namics to determine the k-dependence of the gaps. We show that the optical transmittance, 
conductivity, specific-heat jump, and thermodynamic critical field Hf. are well described by 
a two-band superconductor model with different anisotropics in k-space. The symmetry in 
k-space is determined in order to explain these experiments consistently. 

The optical conductivity for anisotropic s-wave SC is investigated and compared with 
available data for MgB2. A simple angle-dependent generalization of the Mattis-Bardeen 
formula IsgI is used to calculate the optical conductivity. The density of states A^(e) = 
e/Ve^ — is generalized to A^(e) = (Ree/ \/ ^ — A^)^, where the bracket indicates the 
average over the Fermi surface. We employed the following formula at T=0: 

-^H = -/ ciE[Ar(i?)iV(^-E)-(Re^=^).(Re^=^===),,], (5.1) 
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- (Re , :)fc(Re , (5.2) 

We here mention that if the samples are clean and belong to the category of London super- 
conductors, we must use the formulas in the London limit. For a clean superconductor, it 
seems better to use a2s for the two-fluid model|4^. The optical data that we will consider 
here exhibit behaviors explained by the conventional formulas of Mattis and Bardeen. The 
anisotropic order parameters considered in this paper are: 

Aci(k) = A(l + acos(2^)), 
A,2(k) = A(l-6cos2(^^)), 

Aab{k) = A(l + ccos(60)). (5.3) 

Here, 6 and are the angles in the polar coordinate where 6 is the polar angle with respect 
to the c-axis. The parameters a, b and c determine the anisotropy. Ad is a prolate form gap 
for a > and is oblate for a < 0. Ac2 {b > 0) shows the same anisotropy as Ad for a < 0. 
Aab(k) indicates an anisotropy in the a6-plane; the SC gap may possibly be anisotropic in 
the plane since the 2D-like Fermi surface has a hexagonal symmetry. j^O] The integral in 
eq. (j5.ip is evaluated numerically by writing the average over the Fermi surface with elliptic 
functions. For example, for Ad(k) = A(l -|- acos(2^^)), the average over the Fermi surface 
for < a < 1 is given by 



<^Vc^--Ad(k)-^ - W (i' ') ' ~ - " - + 

= ^\/^^(^'^)' (1 + «)A<^, (5.4) 



where P = (u; - (1 - a)A)/(2u;), 7 = sin" V4aAcc;/[(u; - (1 - a)A){uj + (1 + a)A)] and 
F(7, k) is the elliptic integral of the first kind. 

First, we examine a one-band anisotropic model and show that the one-band model is 
insufficient to understand consistently optical and thermodynamic behaviors. In Fig. ITTlthe 
t.a,.„...ion r. at r = is shown as a function of the faouency ..Q We again employ 
the following phenomenological expression for Tq/Tm\1FI l61l|. 

Ik = I (55) 

[Til' + (1 - T'f){ajar.)]^ + [(1 - Tl/'){a,/anW' 
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TABLE II: Several physical quantities obtained by the two-band model with c ~ 0.33 {a band) 
and a ^ 0.33 (vr band). 



Two-band 0.45/0.55 ~ 0.35 ~ 0.82 ~ 0.95 
Exp. 0.45/0.55 0.3 -0.4 0.76- 0.92 0.96 



where cxi and a2 are real and imaginary parts of the optical conductivity, respectively. T/v 
is determined from the expression for the ratio of the power transmitted with a film to 
that transmitted without a film given as = 1/[1 + a,4Zo/{n + Here, d is the 

film thickness, n is the index of refraction of the substrate, and Zq is the impedance of 
free space. We have assigned the following values: d = 10"^ cm, n ^ 3, Zq = 377Q, and 
a„ ~ 8 X lO^r^^^cm^^. Then we obtain T/v — 0.014. The theoretical curves for Ts/Tn 
are shown in Fig. ^2 they have peaks near uj ~ 2Aq. For the oblate, its peak shows an 
increase only twice the normal state value, while the prolate and ab-plane anisotropics show 
more than twofold increases. The experiments show an approximately 2.5-fold increase 2^ 
which supports the prolate or ab-plane anisotropic symmetry. However, the temperature 
dependence of the ratio if^2^/if^2) which increases as , indicates 

that A(k) has an oblate form instead of a prolate form^^l in contrast to Ts/T^. It is also 
difficult to describe the thermodynamic quantities such as the specific-heat jump at T = Tc 
and the thermodynamic critical magnetic field within the single-gap model consistently. 
The specific- heat jump at Tc is given by 

AC(Tc) _ 12 {zr 
7cT. 7C(3) {z^) ' ^'-'^ 

where 7c is the specific-heat coefficient and z is an anisotropic factor of the gap function. 
{z"") is the average of z"' over the Fermi surface. In Fig. ^Jthe specific- heat-jump ratio vs 
anisotropy (a or c) is shown. The experiments indicate that this value is in the range of 
0.76- 0.920,0; the fitting parameters must be a ~ 0.3, —a ~ 0.5 and c ~ 0.3 for the 
prolate, oblate and ab-anisotropic types, respectively. We must assign different values to 
parameters a and c in order to explain the thermodynamic critical magnetic field He- The 
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ratio of Hc{T = 0) to the BCS value is given as 

£^ . ^(.^) . 5.94(.^ (5.7) 

. . ..... j:..0, .... . . . ... 

be small, a ~ 0.07, for the prolate form, and the ab-plane anisotropic and oblate forms 
(a < 0) are ruled out since (2;^) > 1. In Table I, we summarize the status for the single-gap 
anisotropic s-wave model applied to MgB2. As shown here, it is difficult to understand 
the physical behaviors measured using several experimental methods consistently within the 
single-gap model. 

Here, a two-band model with two different anisotropics is investigated. We assume that 
the hybridization between a and tt bands is negligible, and that the optical conductivity is 
given by 

a = W„a„ + ti'vrCrvr, (5.8) 

where and (j,r denote the contributions from a- and vr-bands, respectively. For the case 
of isotropic two gaps, ai must have a shoulder-like structure which appear as an addition of 



two contributions from the two bands, if the magnitudes of two SC gaps are different. 



'he 



experimental data of cxi, however, does not have such a sharp structure (see Fig. IT^ . 
Therefore we must take account of anisotropics for the two-band model. We assume the in- 
plane anisotropy for the two-dimensional-like cr-band, while we assign the three-dimensional 
anisotropy to vr-band where the prolate and oblate forms are examined. 

The transmission Ts/Tn in Fig. IT^ shows that the theoretical curve is in good agreement 
with the experimental curve. The optical conductivity is also described well by the two-band 
model as shown in Fig. We assign the following parameters to the best fit model in Figs. 
El and ^3 the cr-band has ab-plane anisotropy with c ~ or less than 0.33 and the vr-band 
has the prolate form gap (cigar type) with a ^ 0.33. The ratio of the weight of the a-band to 
that of the 7r-band is 0.45/0.55, which agrees with penetration depthjs^ and band structure 
calculations [sl. The ratio of the minimum gap to the maximum gap is 0.35, which is in the 
range of previously reported experimental values. 14 , [s^ Let us mention that the effect of 
cr-band anisotropy is small for the transmission T^/T/y. 

In Fig. El the thermodynamic critical rnagnetic field Hc{T) is shown for the single-band 
and two-band models with available data.p We have simply assumed that the total free 
energy is given by the sum of two contributions from a- and vr-bands: Q = Wa^Q„ + Wt^Q.,,. 
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The experimental behavior is well explained by the two-band anisotropic model using the 
same parameters as those for Tg/T/v and cris/cri„. We show several characteristic values 
obtained from the two-band model in Table II. 

Let us mention here that the two-gap model shows consistency concerning other physical 
quantities. Results of analyses of Hc2 and specific heat usin g th e effective mass approach are 
consistent with those obtained using the two-band model. UQ, 13 been reported that 

the increasing nature of Hc2,ab/Hc2,c with decreasing temperature is explained by the two- 



Fermi surface model. |38|| The specific-heat coefficient 7 in magnetic fields seems consistent 
with that of the multiband superconductor, p, |40| 

VI. SUMMARY 

We have discussed the optical properties in unconventional superconductors. Theoretical 
aspects of the conductivity were discussed in detail from the linear response theory to the 
formula in the London limit. We have presented a new method (R-T method) to measure 
(j{uj) in the far- infrared region from reflectance and transmittance data without the use of 
the Kramers-Kronig transformations. This method provides a method to obtain the far- 
infrared properties more precisely compared to the conventional method. The conductivity 
sum rule is discussed briefly. It has been reported that the sum rule is satisfied for the 
optical conductivity spectra of NbNi-ajCj; that is a typical conventional superconductor. 

We have successfully made a comparison between experiments and theory for the opti- 
cal conductivity of Nd2_2^Ce2.Cu04 in the far-infrared region. We have shown that there 
is a reasonable agreement between the optical conductivity observed by the R-T 

method and theoretical analysis without adjustable parameters except the superconduct- 
ing gap. An estimate of 60~70 cm~^ for the superconducting gap is consistent from both 
the experimental and theoretical aspects. The far-infrared optical conductivity suggests 
that the superconducting gap of electron-doped Nd2_a;Ce2.Cu04 is unconventional one with 
nodes on the Fermi surface. The anisotropic nature of electron-doped superconductors is 
consistent with the recent research performed for the one-band and three-band Hubbard 
models. 0, H B, Q If the superconducting gap is anisotropic for the electron-doped 
superconductors, there is a possibility that both the hole-doped and electron-doped cuprates 
superconductors are governed by a same superconductivity mechanism. 



38 



We have also examined the transmittance, optical conductivity, specific-heat jump and 
thermodynamic critical magnetic field He of MgB2 based on the two-band anisotropic s-wave 
model. This material is described by two order parameters attached to a- and vr-bands, 
respectively, which, moreover, have further anisotropy. We have shown that the two-gap 
model with different anisotropy in k-space can explain the experimental results consistently. 
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